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Abstract 

We consider tiling rectangles of size 4m x 4n by T-shaped tetro- 
minoes. Each tile is assigned a weight that depends on its orientation 
and position on the lattice. For a particular choice of the weights, the 
generating function of tilings is shown to be the evaluation of the mul- 
tivariate Tutte polynomial Zg{Q,v) (known also to physicists as the 
partition function of the Q-state Potts model) on an (m — 1) x (n — 1) 
rectangle G, where the parameter Q and the edge weights v can take 
arbitrary values depending on the tile weights. 

The problem of evaluating the number of ways a given planar domain 
can be covered with a prescribed set of tiles (without leaving any holes or 
overlaps) has a long history in recreational mathematics, enumerative com- 
binatorics, and theoretical physics. 
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Figure 1: Weights of the four types of T-shaped tetrominoes. 
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Figure 2: The domain D to be tiled (here with m = 2, n = 3). Some of the 
vertices are distinguished by black or white circles. 

In this note we consider tiling an M x rectangle of the square lattice 
(henceforth denoted V) with T-shaped tetrominoes, i.e., tiles of size four 
lattice faces in the shape of the letter T. The four different orientations of 
the tiles are shown in Fig. ^ Note that the boundary of each tile comprises 
ten vertices, of which eight are corners. Fig.Elshows the domain T> to be tiled. 
The circles distinguish two subsets of the vertices which we shall henceforth 
refer to as either black (filled circles) or white (empty circles). 

We have the following remarkable 

Theorem 1 (Walkup [IJ). T) is tileable iff M = 4m and N = An. In a 

valid tiling, no tile corner covers a white vertex, whereas black vertices are 
covered by tile corners only. 

Theorem 1 implies that tiles can be distinguished not only by their ori- 
entation (cf. Fig. P) but also by their position relative the the white vertices. 
There are two possible situations, as shown in Fig. depending on which of 
the two cornerless vertices of the tile covers the white vertex. (This applies 
to any of the four orientations of the tiles, although Fig. IHl only illustrates 
the first orientation.) 

Note that the lattice B of black vertices is a tilted square lattice of edge 
length 2^2. This can be divided into two sublattices, -Beven and -Bodd? which 
are both straight square lattices of edge length 4. Our convention is that all 
black vertices on the boundary of V belong to -Bodd- 
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Figure 3: Weights depending on the relative position of a white vertex within 
a tile. 
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Figure 4: Two ways of decorating a tile by an oriented half-arch, depending 
on its position relative to a white vertex. 



We similarly divide the lattice W of white vertices into two sublattices 
Weven and Wodd- Referring again to Fig. we adopt the convention that 
white vertices on the lower and upper boundaries of V belong to Wodd and 
that those on the left and right boundaries belong to VFeven- 

Equivalently, any vertex in i? or belongs to the odd (resp. even) sub- 
lattice whenever its distance from the bottom of V is twice an even (resp. 
odd) integer. 

Finally, we define the generating function 

Fv{{a}Ab})= E rK'n^f (1) 

T£T{V) t£T j=l 

where T{T)) is the set of all T-tetromino tilings of "D, and t is a single tile 
within the tiling T. The weights a^^^^ depend on the orientation o{t) = 
1,2,3,4 of the tile t (cf. Fig. H]) and on the (unique) white vertex w(t) that 
it covers. Moreover, Bj is the number of tiles of position j = 1,2 relative 
to the white vertices (cf. Fig. IH)). Note that the weights bj are not vertex 
dependent. 

Evidently, the tiles covering the white vertices of the boundary B of the 
domain V have their orientation fixed (their long side is parallel to B). It 
follows that those boundary tiles contribute the same a-type weight to any 
tiling. Without loss of generality, we can therefore henceforth set a^l^^ = 1 
for any t such that w(t) G B. 



3 



(1) 



(2) 



(3) 



(4) 



^1 ^1 r^'-^^iV 



(5) (6) (7) (8) 

Figure 5: Arrow conservation at tile junctions. 



It is convenient to illustrate the distinction of Fig. El in another way, by 
decorating the interior of each tile by an oriented half-arch (see Fig. 0]). We 
then have the following 



Lemma 1. In any valid tiling, the orientation of the half-arches is conserved 
across the junctions of the tiles. 

Proof. By rotational and reflectional symmetry, it suffices to consider the 
junction between a tile in the Oi-type position (cf. Fig. [T]) and the neighbour- 
ing tile immediately North-East of it. 

Consider then the situation where the first tile is of the bi type (cf. Fig.|3}. 
The position of the neighbouring tile is constrained by Theorem 1 to have 
the correct location of its white vertex. When the neighbouring tile has its 
long side vertical, there are just two possibilities, shown in Fig. panels 
(1) and (2). The possibility (1) is disallowed since it leaves an untiled hole 
(shown in black), whereas (2) is allowed. When the neighbour tile's long 
side is horizontal, we have the possibilities (3) and (4). But (3) is actually 
disallowed, since, when adding a further tile to cover the lattice face shown 
in black, a white point would be placed in one of the four vertices marked 
by a cross, and neither of these vertices is a valid position for a white point 
according to Theorem 1. In summary, only the possibilities (2) and (4) are 
allowed, and both of those conserve the arrow orientation across the tile 
junction. 

Consider next the situation where the first tile is of the 62 type. A neigh- 
bour tile with a vertical long side leads to the possibilities (5) and (6), but 
(5) is disallowed since it leads to two overlaps (shown hatched). Similarly, 
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when the neighbour's long side is horizontal, one can rule out (7) due to two 
overlaps, whereas (8) is allowed. In summary, only (6) and (8) are allowed 
and compatible with arrow conservation. □ 



Definition 1. Let G = {V,E) be a finite undirected graph with vertex set 
V and edge set E. The multivariate Tutte polynomial of G with parameter 
Q and edge weights v = {ve}e&E is 



where k{A) denotes the number of connected components in the subgraph 



For the relation of the multivariate Tutte polynomial to the more standard 
two-parameter Tutte polynomial Tg^x^u) we refer to Section 2.5 of [2]. In 
the physics literature, Zg{Q,"v) is better known as the partition function of 
the Q-state Potts model in the Fortuin-Kasteleyn formulation [3 . 

In the remainder of the paper we let V = -Seven and E the associated 
natural set of horizontal and vertical edges of length 4. In other words, G is 
an (m — 1) X (n — 1) rectangle G, as referred to in the abstract. Furthermore, 
let w{e) be the natural bijection between edges e & E and white vertices 
w E W \B that are not on the boundary B of V. 

We can now state our main result which is 

Theorem 2. The generating function Fx)({a}, {b}) of T-tetromino tilings of 
the Am x An rectangle T> depends only on the parameters {a} through the 
combinations aiOs and 020,4. For the specific choice of tile weights 



V, 



(2) 



ACE e&A 




for W{e) e We^en \ B 

otherwise 





otherwise 



(3) 



and 



bi = ib^r' = q 



1/4 



(4) 



we have the identity 



g-"/2F^({a},{6}) = ZG(g;v) 



(5) 



with the correspondence between parameters 



Q = {q + q-'f 



Ve = {q + q )xe for e G -E 



(6) 
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Figure 6: A valid tetromino tiling of the domain T> shown in Fig. ^ The 
vertices V entering the definition of the Tutte polynomial (j21) are shown as 
fat solid circles, and the edges of the subset A (1 E are shown as fat lines. 

Proof. By Theorem 1, any white vertex w e is at the junction between 

the long sides of exactly two tiles. If the long side is horizontal, the contri- 
bution of the a-type weight is a'^a^, and if it is vertical the contribution is 
0^04 . This proves the first part of the theorem. (Above we have already 
explained that one may set all = 1 when w & W (iB.) 

By Lemma 1 the oriented half-arches form continuous curves with a con- 
sistent orientation (see Fig. [Ul for an example). Moreover, by the consistency 
of the position of the white vertices (cf. Figs. I2H21), these curves cannot end 
on B. Therefore, the half-arches form a set of oriented cycles. 

Let C be the set of all possible configurations of cycles arising from valid 
tilings of V, but disregarding the orientations of the cycles. The summation 
in the generating function ((T)) can then be written ^^-(j,) = J2c Sr(X')|C' 
where the last sum is over cycle orientations only. A close inspection of 
panels (2), (4), (6) and (8) in Fig. El reveals that the orientation of a given 
oriented cycle cq can be changed independently of all other oriented cycles in 
the tiling T G T(V) by shifting the tiles traversed by cq by one lattice unit to 
the left, right, up or down, but without moving any tiles not traversed by Cq. 
In particular, this transformation does not change the orientation of any tile 
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(cf. Fig. P), and so does not change the a- type weights. We have therefore 



M{ah{b})=j:u<s E n^f • (7) 

C T(V)\Cj=l 

[The remainder of the proof parallels a construction used in Q to show 
the equivalence between the Potts and six-vertex model partition functions.] 

The partial summation X]r(©)|c 11^=1 appearing in ((7j) can be carried 
out by noting that the 6-type weights (jH) simply amount to weighting each 
turn of Co through an angle ±7r/2 by the factor g^^/"^. Since the complete 
cycle turns a total angle of ±2ti depending on its orientation, the product of 
the 6-type weights associated with the cycle cq, summed over its two possible 
orientations, yields the weight q + q^^ = Q^l"^ . Denoting by ^(C) the number 
of cycles in C, we have therefore found that 

Y.i\hf=Cfl\ (8) 

r(i5)|cj=i 

It now remains to sum over the un-oriented cycles C . There is an easy 
bijection between C and the edge subsets ACE appearing in (j21). Namely, 
an edge e E E belongs to A iff it is not cut by any cycle in C (see Fig. E}. 
Conversely, given A, the cycles are constructed so that they cut no edge in 
A, cut all edges in E\A, and are reflected off the boundary B in the vertices 



W nB. Under this construction, the a-type weights Q simply give Yl 
Collecting the pieces this amounts to 

ACE eeA ^ 

Using finally the Euler relation, and remarking that \V\ = mn, one ob- 
tains ©. □ 

Let us discuss a couple of special cases of Theorem 2. First, when there 
are no a-type weights (i.e., setting a^l^^ = 1 for any o{t) and w{t)), the Tutte 
polynomial has Vg = q + q~^ for any e E E. The tiling entropy defined by 

SG{q)= lim (Ft,)^ (10) 

m,n— >oo 
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is then related to the one found for the Tutte polynomial {alias Potts model) 
by Baxter j3] using the Bethe Ansatz method. For < Q < 4 one sets q = e*^ 
with fi e (0, 7r/2) and finds 



log Sciq) = 



dt 



sinh [(tt — yu)t] tanh(yut) 



(11) 



t sinh(7rt) 



whereas for Q > 4 one sets q 



e and finds 




—nX 



tanh(riA) 



logS'G(g) 



n 



(12) 



n=l 



Finally, for Q = 4 one has 



o rn-rr(5/4)y 



(13) 



An even more special case of the identity ^ arises when 6i = 62 = 1 as 
well. F-p is then the unweighted sum over all T-tetromino tilings of T>. In 
terms of the standard Tutte polynomial Tg{x, y) one then has F-p = 2Tg{S, 3), 
as first proved by Korn and Pak 

References 

[1] D.W. Walkup, Amer. Math. Monthly 72, 986-988 (1965). 

[2] A.D. Sokal, in Surveys in Combinatorics, 2005, edited by Brid- 
get S. Webb (Cambridge University Press, 2005), pp. 173-226; 
|niath.CO/0503607. 

[3] CM. Fortuin and P.W. Kasteleyn, Physica 57, 536-564 (1972). 

[4] R.J. Baxter, S.B. Kelland and F.Y. Wu, J. Phys. A 9, 397-406 (1976). 

[5] R.J. Baxter, J. Phys. C: Solid State Phys. 6 L445-L448 (1973). 

[6] M. Korn and I. Pak, Theor. Comp. Science 319, 3-27 (2004). 



8 



